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ABSTRACT 

We present first results of a Chandra X-ray observation of the rare oxygen- 
type Wolf-Rayet star WR 142 (= Sand 5 = St 3) harbored in the young, heavily- 
obscured cluster Berkeley 87. Oxygen type WO stars are thought to be the most 
evolved of the WRs and progenitors of supernovae or gamma ray bursts. As part 
of an X-ray survey of supposedly single Wolf-Rayet stars, we observed WR 142 
and the surrounding Berkeley 87 region with Chandra ACIS-I. We detect WR 142 
as a faint, yet extremely hard X-ray source. Due to weak emission, its nature as 
a thermal or nonthermal emitter is unclear and thus we discuss several emission 
mechanisms. Additionally, we report seven detections and eight non-detections 
by Chandra of massive OB stars in Berkeley 87, two of which are bright yet soft 
X-ray sources whose spectra provide a dramatic contrast to the hard emission 
from WR 142. 

Subject headings: stars: individual (WR 16, WR 142 = Sand 5 = St 3, BD+36 
4032, HD 229059, V439 Cyg) — stars: Wolf-Rayet — open clusters and associa- 
tions: individual (Berkeley 87) — X-rays: stars 
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1. Introduction 

Wolf-Rayet (WR) stars are massive, highly-evolved stars nearing the end of their lives as 
supernovae (SN) or as collapsing objects emitting gamma-ray bursts (GRB, e.g. MacFadyen 
& Woosley 1999, Postnov & Cherepashchuk 2001, Georgy et al. 2009). WR stars undergo 
rapid mass loss through strong winds, with initial masses of >25M0 (Crowther 2007). The 
classification of WR stars is determined spectroscopically in the optical and is divided among 
the nitrogen-rich WNs, carbon-rich WCs, and oxygen-rich WOs. For the most part, single 
WR stars are thought to follow the evolutionary path: O (LBV/RSG) WN WC 
^ WO ^ SN [or GRB] (Conti et al. 1983, Crowther 2007). The initial mass of the O star 
determines whether it passes through an intermediate luminous blue variable (LBV) or red 
supergiant (RSG) phase (Crowther 2007). 

Sanduleak (1971) noticed a class of stars that did not have planetary nebulae (FN) but 
displayed a WR-like spectrum with strong O VI doublet emission similar to those found in 
the central stars of FN named the O VI sequence by Smith and Aller (1969). These stars 
were suggested to be a separate WO sequence of Wolf-Rayet stars (rather than an extension 
of the WC sequence) by Barlow & Hummer (1982). WO stars are thought to be in the 
late helium-burning stage, or possibly the carbon-burning stage, where the enhanced oxygen 
abundance compared to less evolved WR stars is revealed by mass loss stripping (Barlow & 
Hummer 1982). Of the 298 galactic WR stars in the appended Vllth catalogue of galactic 
Wolf-Rayet stars (van der Hucht 2001, 2006), only four are of the rare WO spectral type, 
including the star WR 142, also known as Sand 5 (Sanduleak 1971) and St 3 (Stephenson 
1966). Crowther et al. (1998) developed a new classification scheme using primary and 
secondary oxygen line ratios that confirmed previous classifications of WR 142 as a member 
of the subclass W02 (Barlow & Hummer 1982, Kingsburgh et al. 1995). 

Despite new discoveries of WR stars in the galactic plane, including a WO type as the 
exciting star of the planetary nebula Th 2-A (Weidmann et al. 2008), little is yet understood 
about the mechanisms that drive high-energy processes such as X-ray emission in single WR 
stars, particularly so for the rare WO stars. Although WR stars have much stronger winds 
and are more evolved chemically, the line-driven instabilities that are thought to give rise to 
soft X-rays (kT < 1 keV) from shocks in O star winds may also be present in WR winds 
(Gayley & Owocki 1995). If that is the case, then WR stars may also be capable of producing 
soft X-rays from radiative wind shocks (Baum et al. 1992). 

Few single WR stars had been studied with high sensitivity in X-rays until recently. 
Thus far, several apparently single WN stars have been observed to emit X-rays (Skinner et 
al. 2002, Ignace et al. 2003, Oskinova 2005). An ongoing Chandra a,nd XMM-Newton survey 
has recently detected the apparently single WN stars WR 2, WR 18, WR 79a, and WR 134 
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with X-ray luminosities log Lx ~ 32.2 - 32.7 erg s ^ (Skinner et al. 2010b), comparable to 
some WR+OB binaries like 7^ Vel (WC8+07) with log Lx = 32.9 erg s"^ (Skinner et al. 
2001) or WR 147 (WN8+0B) with log Lx = 32.83 erg s'^ (Skinner et al. 2007; see also Sec. 
4.4). Only upper limits of log Lx < 29.82 - 30.97 erg s~^ exist from observations of single 
WC stars (Oskinova et al. 2003, Skinner et al. 2006). The WC sequence is thus either faint 
in the X-rays or possibly X-ray quiet. WR 142 just recently became the first WO star to be 
detected in the X-rays using XMM-Newton (Oskinova et al. 2009). 

WR 142 resides in the open cluster Berkeley 87. Table 1 summarizes the general prop- 
erties of WR 142. Berkeley 87 lies in a heavily obscured region of the Milky Way in Cygnus. 
Initial cluster age estimates of ~ 1 - 2 Myr (Turner & Forbes 1982, hereafter TF82) have now 
been revised upward. A study restricted to three of the highest mass cluster stars suggests 
a slightly older age of ~ 3 Myr (Massey et al. 2001). A more recent study (Turner et al. 
2010) gives ~ 5 Myr. Berkeley 87 is an interesting cluster containing 105 cluster members 
identified by an optical study (TF82), B-supergiants including HD 229059, an 08.5-09 star 
BD-l-36 4032. a possible luminous blue variable Be star V439 Cyg, and the pulsating M3- 
supergiant BC Cyg. Additionally, OH masers and compact HII regions trace massive star 
formation 3' - 9' north of WR 142 (Argon et al. 2000, Matthews et al. 1973). At a distance 
of 1230 ± 40 pc (Turner et al. 2006), the proximity of Berkeley 87 offers ample opportunity 
to learn about the properties of rare objects such as WR 142 as well as the surrounding 
cluster. 

As part of an X-ray survey aimed at determining if single Wolf-Rayet stars that are 
not known to be in binary systems are X-ray emitters (specifically which spectral types), we 
observed WR 142 with Chandra. We are unaware of any evidence pointing to a companion of 
WR 142. WR 142 was chosen over other WO stars because of its relatively nearby distance, 
low Av compared to other WR stars, and interesting surrounding cluster. Its astonishing 
supersonic wind (voo = 5500 km s~^, Kingsburgh et al. 1995) and an unusual optical detection 
of diffuse emission from the C IV doublet (Polcaro et al. 1991) that is a signature of shocked 
gas make WR 142 an even more interesting target as one of the few oxygen-type WR stars. 

The primary objective of the study presented here is to determine the X-ray properties of 
WR 142. X-ray emission from WR 142 was previously observed by XMM-Newton (Oskinova 
et al. 2009), but Chandra has several advantages over XMM-Newton, including improved 
angular resolution (by a factor of ~ 4 - 5), sharp on-axis point-spread function (PSF) to 
provide checks on source extent, longer uninterrupted exposure that provides a continuous 
light curve and more stringent test for variability, and lower instrumental background. Such 
benefits allow for reliable source identification, especially when compared with optical po- 
sitions; negligible background subtraction for lightcurves and spectra; and searches for the 
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presence of diffuse shock emission around WR 142. We analyze tlie faint X-ray spectrum of 
WR 142 and discuss possible X-ray emission mechanisms. We also report Chandra detections 
of seven other massive OB stars in Berkeley 87 and compare the much softer spectra from 
two X-ray bright OB stars to the hard spectrum of WR 142. In this paper, we focus on WR 
142 and massive OB stars in Berkeley 87. A first overview of the Chandra results for the 
cluster as a whole was presented by Skinner et al. (2010a) and the cluster will be discussed 
in more detail in a future paper. 



2. Chandra Observations and Data Reduction 

The Chandra observation (ObsID 9914) centered near WR 142 began on 2009 February 
2 at 02:41:24 TT and ended on February 2 at 22:52:55 TT. The nominal pointing coordinates 
were located 28" NE of WR 142 at R.A. = 20/i 21m 46.42s, Decl. = +37° 22' 44.9" (J2000.0). 
The observation was obtained with the ACIS-I (Advanced CCD Imaging Spectrometer) 
imaging array in faint timed-event mode with 3.2 s frame times. The exposure live time 
was 70148 s. ACIS-I consists of four front-illuminated 1024 x 1024 pixel CCDs with pixel 
size ~ 0.492" and has a combined field of view (FOV) of ~ 16.9' x 16.9'. More information 
on Chandra and its instrumentation can be found in the Chandra Proposer's Observatory 
Guide (POG)@. 

The Chandra X-ray Center (CXC) provided a Level 2 events file, which was analyzed 
with CIAO version 4.1.20 using standard science threads. To determine the source centroid 
position, the net counts, and a 3cr extraction region of the target WR 142, we used the CIAO 
wavdetect tool. We ran wavdetect on full-resolution images using events in the 0.3 - 8 keV 
range, chosen to reduce background. The wavdetect threshold was set at sigthresh = 1 x 
10~^. Scale sizes of 1, 2, 4, 8, and 16 were used. The 3a extraction region determined for WR 
142 by wavdetect was an ellipse with an area of ~ 11.25 pixels^ (Figure 1) and was used for 
spectral and timing analysis. The CIAO task srcextent, which utilized point- spread-function 
(PSF) information, was also used to determine the source size of WR 142 and determine if 
its emission is point-like or extended. 

Within the 0.3 - 8 keV energy range, light curves were extracted using CIAO dmextract 
and source variability was tested by two different methods. Using unbinned photon arrival 
times, we applied the Kolmogorov-Smirnov (KS) test (Press et al. 1992). For comparison. 



^See [http : / /asc. harvard.edu/proposer /POG| 

^Further information on Chandra Interactive Analysis of Observations (CIAO) software can be found at 
http://asc.harvard.edu/ciaol 
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the Gregory-Loredo algorithm (Gregory & Loredo 1992, 1996) was apphed by using CIAO 
(?/war?@ tool. 

CIAO specextract was used to extract source and background spectra. The tool specex- 
tract also created response matrix files (RMFs) and auxiliary response files (ARFs). Spectral 
analysis utilized the HEASOFT Xana(if^ software package including XSPEC vers. 12.5.0. 
The background region was chosen to be a source-free annulus centered on WR 142 of r = 4 
to r = 16 pixels. From this region, we estimate < 1 background count (0.3 - 8 keV) inside the 
3a elliptical source region for WR 142. The background inside the source region is negligible 
relative to the number of source counts. 



3. Results 

3.1. Imaging and Source Identification 

WR 142 is detected by Chandra as an X-ray source at J202144.35+372230.7 (Figure 1). 
Table 2 summarizes some properties of the X-ray emission from WR 142. These coordinates 
are the source centroid as determined by wavdetect using an energy filtered image (0.3 - 8 
keV). Without further astrometric calibration, the Chandra ACIS-I absolute position has an 
accuracy of ~ 0.4" (68% limit j^. The X-ray position of WR 142 is in exact agreement with 
the 2MASS position of J202144.35+372230.7. Additionally, the C/ianc/ra position has offsets 
from optical counterparts of only 0.20" from the USNO-Bl.O source at J202144.35+372230.9 
and 0.23" from the HST GSC2.3.2 source at J202144.36+372230.5. All offsets are well within 
the positional accuracy of Chandra. 

Results from the CIAO tool srcextent indicate WR 142 is not extended. Use of a PSF file 
specifically tailored for our observation with the WR 142 spectrum and the position of WR 
142 on the ACIS-I detector yields a source size of only 0.52" (0.41-0.64", 90% confidence), 
equivalent to about one ACIS pixel, and a PSF size of 0.85" (0.66-1.03", 90% confidence). 
Thus, the X-ray emission of WR 142 is not extended at the 90% confidence level. 



^http://cxc. harvard.edu/ciao/ahelp/glvary.litml. 

^The XANADU X-ray analysis software package is developed and maintained by 
NASA's High Energy Astrophysics Science Archive Research Center (HE AS ARC). See 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xanadu.html for further information. 



^http://cxc.harvard.edu/cal/ASPECT/celmon/. 
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3.2. Spectral Analysis 

The X-rays from WR 142 are weak but extremely hard, with a mean photon energy of 
4.31 keV (Table 2). Surprisingly, WR 142 does not exhibit any significant emission below 
2 keV and there is no evidence of a soft component (Figures 2, 3). We find only one event 
below 2 keV, with an energy of 1.89 keV. This is in contrast to a faint soft component at 
E < 1 keV in the XMM-Newton EPIC MOS spectrum of WR 142 reported by Oskinova et 
al. (2009). Assuming that no variabihty occurred in the WR 142 source spectrum or Nh 
between the two observations, the most likely explanation for this difference would be the 
much lower effective area of Chandra ACIS-I at low energies E < 1 keV. 

We tried to fit the spectrum of WR 142 (Figure 3) with several different models: a 
IT (one temperature) APEC optically thin plasma model and a power-law model. The 
APEC model requires a very high temperature and thus a thermal bremsstrahlung model is 
similar. Because of low counts, the spectral fit results were somewhat sensitive to binning 
strategy. Representative spectral fits are shown in Table 3, based on fits of unbinned spectra. 
Between the different models, there is little variation in the quality of fit using C-statistics 
(Cash 1979). However, determining whether a particular model is statistically acceptable 
using C-statistics is not as straightforward as with statistics (Heinrich 2003). 

Allowing the abundances to vary from solar displays no significant improvement to the 
APEC fits, nor does the use of representative WO abundances (Table 3). This is largely 
because we do not detect any emission lines in the spectrum from elements such as O which 
are highly nonsolar in WO stars. Additionally, all of the O and C lines that lie in the Chandra 
bandpass are below 1 keV where we do not detect any signal. 

The data show a slight preference for a model that includes a gaussian line for an iron 
line in the ~ 6.4 - 6.67 keV energy range (Table 3). The counts from WR 142 are not 
sufficient to distinguish between the possible iron lines. The event energies closest to the 
possible iron emission, 6.24, 6.40, 6.60, and 6.63 keV, show two events within 0.1 keV of the 
iron K line complex at 6.67 keV (Fe XXV He-like triplet). If a weak Fe line is present, it is 
likely to be the 6.67 keV iron K line, which has been detected in other WR stars. 

We have also experimented with fits of binned spectra using statistics in order to 
estimate confidence bounds on important fit parameters (Fx, Nh, kT, Fpow)- Given the 
limited number of counts (46), binning at desired levels of >10 counts per bin for fitting 
with statistics was not feasible. We thus consider spectra binned to a minimum of 6 
counts per bin. Results indicate that if the plasma is thermal, as the APEC models assume, 
then the temperature kT must be very high, with a lower 90% confidence bound of kT > 
3.6 keV (Table 3). The Chandra data do not place a useful upper bound on kT. 
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In determining the X-ray flux, the unbinned spectral fits display fairly consistent results. 
The fit with the lowest C-statistic shown in Table 3 gives an observed (absorbed) X-ray flux 
of Fx = 1.8 (0.8 - 2.5) X 10~^^ erg s~^ cm~^ (0.3 - 8 keV, la confidence intervals). The error 
bounds reflect the range of fluxes measured for the different models in Table 3. The best fit 
models give X-ray luminosities of log Lx = 30.72 - 30.82 (ergs s~^), which are derived from 
unabsorbed fluxes (0.3 - 8 keV) resulting from different unbinned spectral fits in Table 3. 

All indications are that the absorption towards WR 142 is high, with Nh around 4 - 
5 X 10^^ cm^^. This is about 3-4 times higher than the absorption estimated from the 
observed optical extinction of 1.3 x 10^^ cm~^ obtained using Ay ~ 6.1 (van der Hucht 2001, 
Gorenstein 1975 conversion). Spectral fits of WR 142 with Nh fixed at the value determined 
from optical extinction give much larger fit residuals than the fits with larger absorption 
given in Table 3. 

3.3. Timing and Vciriability Analysis 

The source shows no obvious large-scale variations nor flares in its X-ray lightcurve 
(Figure 4). As for short term variability, KS statistics and the Gregory-Loredo algorithm 
give a probability of constant count rate of 0.77 and 0.40, respectively. Thus, neither test 
imphes variabihty at high confidence levels (>90%) on timescales of hours during the <1 
day Chandra observation. 

For purposes of comparing the observed X-ray flux by Chandra and the previously 
pubhshed XMM-Newton X-ray flux of Fx = 4 ± 2 x 10"^^ erg cm-^ s"^ from 0.25 - 12 keV 
(Oskinova et al. 2009), we measure Fx = 3 ± 1 x lO"^'' erg cm-^ s"^ (0.25 - 12 keV) using 
an average value of the fluxes from the extended energy range from the three models shown 
in Table 3. Thus, Chandra and XMM-Newton X-ray fluxes agree to within the uncertainties. 
Based on the above flux comparison, any variability over the approximate one year interval 
between the Chandra and XMM-Newton observations was a factor of ~ 2 or less. 



3.4. Summary of WR 142 X-ray Properties 

Based on the Chandra data discussed above, the X-ray properties of WR 142 can be 
summarized as follows: (i) source structure is consistent with a point source at Chandra^s 
angular resolution, (ii) no significant variability is detected over a timescale oi 1 day, (iii) 
strong X-ray absorption below 2 keV, with an absorption column density (Nh) that is a 
factor of ~ 3 - 4 greater than expected from published optical estimates, (iv) no evidence 
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for a soft X-ray component below 2 keV, (v) a spectrum dominated by hard X-rays above 
2 keV that can be fit equally well using a very high temperature thermal plasma model or 
a power-law model, and (vi) a relatively low X-ray luminosity log Lx(0.3 - 8 keV) = 30.7 - 
30.8 ergs s~^ (at d = 1.23 kpc), which is ~ 10 times less than found for two X-ray bright 
OB stars in Berkeley 87 (Section 4). We discuss these properties in more detail below. 

4. OB Stars in Berkeley 87 

In this section, wc comment briefly on OB stars in Berkeley 87 contained in the Chandra 
FOV. We report Chandra detections of seven OB stars and identify eight undetected B-type 
stars (Table 2) that are likely cluster members (TF82). Figure 1 shows the relative positions 
of the X-ray bright stars and the noteworthy non-detection of the LBV candidate V439 Cyg. 

4.1. BD+36 4032 

This star was classified as 09V by TF82 (Berkeley 87-25) and was the only object 
explicitly identified as an 0-type star in their catalog of 105 likely members of Berkeley 87. 
The spectral type was later refined to 08.5 111 by Massey ct al. (2001), who also estimated 
the mass to be = 39 Mq. BD+36 4032 hes ^ 3' northwest of WR 142. 

X-ray characteristics are listed in Table 2. The X-ray spectrum is soft (Figure 3) with 
almost all of the emission emerging below 2 kcV, which is dramatically different than the 
hard emission from WR 142. Reasonably good fits can be obtained using a solar abundance 
IT APEC model with a cool plasma component at kT 0.6 keV (Table 4). The IT APEC 
fit can be improved slightly by adding a second temperature component between 1-2 keV, 
but almost all of the emission measure still resides in the cool component. 

The best-fit Nh given in Table 4 equates to Ay = 5.1 [4.4 - 6.5] mag (Gorenstein 1975). 
This is in good agreement with the optically-determined value Ay = 4.8 mag using E(B — 
V) = 1.60 and Ay = 3E(B - V) (TF82). Thus, the excess X-ray absorption detected in WR 
142 is not present toward this O star. 

4.2. HD 229059 

This B supergiant has V = 8.71 mag and was the brightest member of Berkeley 87 listed 
by TF82 (Berkeley 87-3). They assigned a spectral type B2 labe but Massey et al. (2001) 
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classified it as Bl la. 

The Chandra spectrum of HD 229059 is nearly identical to that of BD+36 4032 (Figure 
3). One noticeable difference is that HD 229059 shows stronger emission in the 1.8 - 1.9 keV 
range, which is likely due to stronger hue emission from the Si XIII triplet (1.84 - 1.87 keV). 

Spectral fits of HD 229059 give nearly identical results to BD+36 4032. Table 4 summa- 
rizes the best-fit IT APEC model. No significant fit improvement was obtained by adding 
a second temperature component. The Nh from Table 4 gives Ay — 5.2 [4.5 - 5.6] mag, in 
excellent agreement with the optically-determined value Ay = 5.1 mag (TF82). Thus, we 
find no excess X-ray absorption toward HD 229059. 

4.3. Berkeley 87-4 

This BO. 2 HI star lies ^ 5.1' west of WR 142 and is object number 4 in the hst of hkely 
cluster members compiled by TF82. It was included in the study of Massey et al. (2001) 
who obtained a mass estimate M* = 23 Mq. Berkeley 87-4 was the third brightest Chandra 
X-ray detection of the OB stars in Berkeley 87 (Table 2). There are sufficient counts for a 
crude spectral fit, and a IT APEC model gives Nh and kT values that are similar to the 
brighter OB stars BD+36 4032 and HD 229059 but a lower X-ray luminosity (Table 4). The 
X-ray derived Nh agrees to within the model uncertainties with that expected from optical 
reddening (TF82). 

4.4. V439 Cyg 

This variable star is also known as MWC 1015 and is object number 15 in the list of likely 
cluster members compiled by TF82. V439 Cyg has anomalous colors and is displaced well 
to the redward side of the zero-age main-sequence in the cluster color-magnitude diagram 
(TF82). The star's spectrum is peculiar and evidently variable, with most recent studies 
assigning an early B spectral type of BOe (Polcaro & Norci 1998) or B[e] (Massey et al. 
2001). V439 Cyg shows some similarities to a LBV (Polcaro & Norci 1998). 

V439 Cyg lies ^ 3.2' NW of WR 142 and is a bright near-lR source with 2MASS data 
giving Kg = 7.59 mag and position 2MASS J202133. 57+372451. 57. There is no X-ray source 
at or near this position in our Chandra image. A circular region of radius 1.96" (90% encircled 
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energy at E = 2 ke\0 centered on the V439 Cyg 2MASS position encloses only one photon 
in the ACIS-I image (0.3 - 8 keV range). 

To compute an X-ray upper limit, we assume a conservative 4-count ACIS-I detec- 
tion threshold and an intrinsic thermal stellar spectrum similar to the two OB stars dis- 
cussed above. Using Nh = 1.0 x 10^^ cm~^ corresponding to E(B— V) = 1.53 (TF82) and 
a Raymond-Smith solar-abundance thermal plasma model with kT = 0.6 keV, the Chandra 
PIMmS simulator gives an unabsorbed flux upper limit Fx(0.3 - 8 keV) < 3.9 x 10 ergs 
cm~^ s~^. If V439 Cyg is a cluster member, then d = 1.23 kpc gives log Lx < 29.85 (ergs 
s-i). 

Given the unusual nature of V439 Cyg, further consideration of its X-ray non-detection 
is warranted. In this regard, it is important to note that Chandra detected only the most 
luminous B stars in the cluster having absolute magnitudes my < —2.1 (Sec. 4.5). Several 
less- luminous early B stars were undetected (Table 2 notes). This is likely a consequence of 
the known Lx oc Lboi trend for O and early B-type stars (Berghofer et al. 1997). The absolute 
magnitude of V439 Cyg is not well-known due to uncertain reddening, but it appears to lie at 
or near the my cutoff for our Chandra detections (Figure 7 of TF82). Thus, V439 Cyg may 
have insufficient Lboi to produce an X-ray luminosity above the Chandra detection threshold. 

Excess absorption toward V439 Cyg may have also partially contributed to its non- 
detection. The anomalous colors of V439 Cyg imply excess reddening, and it may be sur- 
rounded by a shell of material that was ejected during a previous mass loss phase (Polcaro 
& Norci 1998). A gas-rich shell would increase the line-of-sight photoelectric absorption to 
lower energy X-rays (E < 1 keV). 

If V439 Cyg is indeed a LBV (Polcaro & Norci 1998), it could be evolving into the WN 
phase. For stars of initial masses ~ 40 - 75 Mq, the evolutionary sequence is (Crowther 
2007): O LBV WNH-poor- Some putatively single WN7-9 stars (also known as WNL 
stars) have so far gone undetected in X-rays (Figure 5). But a few recent WNL detections 
have been reported, including an unambiguous detection of the apparently single WN9ha 
star WR79a (Skinner et al. 2010b), a possible faint detection of WR 16 (WN8h) in a recent 
XMM-Newton observation (Figure 5), and a new Chandra detection of the WN8 component 
of the close WR+OB binary system WR 147 (Zhekov & Park 2010). Although the nature 
and evolutionary status of V439 Cyg remain somewhat of a mystery, its X-ray non-detection 



^^http://cxc. harvard.edu/cal/Acis/CaLprods/psf/eer_on. html. 

^^For information on PIMMS (Portable Interactive Multi-Mission Simulator) see 
http: / / cxc.harvard.edu/ciao/ahelp/pimms.html 
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would not be in discord with an evolved star that is approaching the WNL phase, given that 
a few other WNL stars have also gone undetected in X-rays. 

4.5. Other B Stars in Berkeley 87 

In addition to the OB stars discussed above, Chandra provides faint X-ray detections of 
four other B stars in Berkeley 87. These are Berkeley 87 no. 13 (B0.5 III), no. 24 (Bl lb), 
no. 26 (BO. 5 lab), and no. 32 (BO. 5 III). Their basic properties are summarized in Table 2, 
but they lack sufficient counts for spectral analysis. We also report non-detections of seven 
other B stars in Berkeley 87 (Table 2 notes). 

The Chandra results show a noticeable trend in that all of the B star detections have 
high luminosities (my < —2.1; TF82 photometry). Most, but not all, of the undetected B 
stars are less luminous (my > —2.1). However, there are a few interesting exceptions. The 
Bl V star Berkeley 87-31 (my = —2.7; TF82) lies nearly on-axis at an offset of 18.7" from 
WR 142 but was not detected. However, it does lie near an ACIS-I CCD gap. Two other 
luminous but undetected B stars, no. 9 (BO. 5 V) and no. 38 (B2 HI), lie further off-axis 
(3.9' and 5.0', respectively). 

5. Discussion 

5.1. Excess X-ray Absorption Toward WR 142 

Chandra spectral fits of WR 142 (Table 3) give an absorption column density of Nh ~ 
4 - 5 X 10^^ cm~^. In contrast, optical studies yield Nh ~ 1 x 10^^ cm~^ (TF82, Barlow 
& Hummer 1982, van der Hucht 2001, Gorenstein 1975 conversions) for Berkeley 87 cluster 
members. Our X-ray analysis of the cluster members BD-l-36 4032 and HD 229059 (Sec. 4) 
also yields Nh ~ 1 x 10^^ cm~^. Thus, the X-ray spectra clearly show that excess absorption 
is present toward WR 142 that is not seen toward two other massive cluster members. This 
could arise either in the dense WO wind or perhaps in cold circumstellar gas. 

Excess X-ray absorption seen as a disagreement between absorption based on optical 
extinction and the X-ray fit Nh has been seen in other WR stars as well, such as 7 VeP 
(Skinner et al. 2001) and recently detected WN stars (Skinner et al. 2010b). 

Inhomogeneous extinction, although present, does not account for the discrepancy in X- 
ray and optical Nh values for WR 142. The spread in extinction values for cluster members is 
small, from Nh = (0.8 - 1.3) x 10^^ cm"^ (TF82, Gorenstein 1975 conversion). Additionally, 
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the smallest foreground extinction is near the south and WR 142 is ~ 2.6' southeast of the 
cluster center (TF82), indicating the excess Nh is due to local absorption associated with 
WR 142. 

If the excess absorption is due to the wind, we can estimate the radius at which the 
X-rays emerge using a wind optical depth calculation. From the observed Chandra WR 142 
spectrum, the X-rays are entirely absorbed below ~ 2 keV. Using generic WO abundances 
(van der Hucht et al. 1986) and assuming the mass- loss parameters in Table 1, the radius 
of optical depth unity at E = 2 keV is R(r = 1, 2 keV) 1500 R©. Clumping in the WR 
winds may reduce the mass loss rate by a factor of 2 - 4 (Crowther 2007) and would reduce 
R(t = 1,2 keV) by the same factor. R(r = 1) is well outside the wind acceleration zone and 
the wind will already have reached terminal speed, assuming a standard /3 = 0.5 - 1.0 wind 
velocity law. It should be kept in mind that this is a minimum radius for 2 keV X-rays to 
escape. 

Another possibility regarding excess absorption is the presence of cold circumstellar 
gas. Nebulosity around WR 142 has been detected (Miller & Chu 1993), though a ring 
nebula is not present. Polcaro et al. (1997) note excess reddening surrounding WR 142 
of at least Nh = 2.7 x 10^^ cm'^ (Gorenstein 1975 conversion) suggesting dense material, 
possibly a nebula from the stellar wind. Alternatively, Lozinskaya (1991) detected a distant 
IR shell (~ 23 pc from WR 142) whose dynamical age suggests that it formed before, and 
therefore independently, of the WO wind. Though it is apparent that some circumstellar 
material is present, our X-ray spectral fits have not distinguished between wind absorption 
and absorption by circumstellar gas far from the star. Both the wind and local circumstellar 
material could contribute to the excess X-ray absorption. 



5.2. Is WR 142 a Binary? 

There is no evidence from ACIS images (or other observations) that WR 142 is a binary. 
For an on-axis point source, the Chandra HRMA and ACIS PSF encircles 70% of the source 
energy for an angular source radius of ~ 1.2" at 2 keV and 1.3" at 6 keV0. If two X-ray 
sources are present, their separation must be <1". 

The WR 142 X-ray luminosity from this Chandra observation is much less than that of 
known WR+0 binaries, which are generally very luminous in the X-rays. Due to plasma 
formed in colliding wind shocks, these systems can have luminosities approaching Lx ~ 10^^ 



http://cxc.harvard.edu/cal/Acis/CaLprods/psf/eer_on. html. 
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ergs s-^ such as 7^ Vel (WC8+07.5III; Skinner et al. 2001, Schild et al. 2004). WR 142 is 
at least two orders of magnitude less luminous. 

However, the possibility of an unseen companion remains. But, a compact companion 
seems to be unhkely. Theoretical models of the WR winds accreting onto neutron stars 
predict high luminosity log Lx > 34 (ergs s"^) (HD 50896 = EZ CMa, Stevens & Willis 
1988). Additionally, a neutron star would be required to spin fast enough to throw off most 
of the accreting material to remain undetected (propeller effect, e.g. Lipunov 1982). If an 
optically faint lower mass normal stellar companion (rather than a compact object) lies close 
to WR 142, it could easily escape detection. 

Currently, the binary frequency for galactic WR stars is ~40% (van der Hucht 2001). 
Additional searches for binarity are needed for WR 142. These could include high resolution 
infrared imaging or radio interferometry, though they would only be capable of detecting a 
companion down to separations of a few tenths of an arcsecond. More promising approaches 
would be to search for periodic photometric or spectroscopic variability that could signal a 
closer spectroscopic companion. Spectroscopically, an X-ray emitting companion orbiting in 
the WR wind could reveal itself through phase variability of the WR star ultraviolet hues 
(Hatchett & McCray 1977). 

5.3. Thermal versus Nonthermal Emission 

Low counts in the WR 142 ACIS-I spectrum result in nearly identical fits by either ther- 
mal APEC models or power-law models. As such, the WR 142 spectrum cannot sufficiently 
distinguish between a very high temperature thermal plasma and nonthermal processes e.g. 
inverse Coulomb scattering due to relativistic electrons accelerated in wind shocks (Chen & 
White 1991a). In general, other WR stars with higher signal-to-noise spectra show detec- 
tions of emission lines indicative of thermal emission. While the goodness-of-fit for spectral 
models slightly improves with the addition of the Fe K line for WR 142, there are not enough 
line counts to qualify as a definite line detection. There is no obvious evidence of other ther- 
mal emission lines, such as the S line at 2.46 keV. Theoretical models do exist that attempt 
to account for the hard X-ray emission from some early type stars by nonthermal processes 
(Chen & White 1991a). However, if WR 142 were a nonthermal X-ray source, it would be 
the first putatively single WR star in that class. The possibility of nonthermal emission 
has been previously mentioned for other WR stars such as WR 110 but the evidence for 
nonthermal emission was not conclusive (Skinner et al. 2002). 
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5.4. Thermal Emission Mechanisms 

Thermal X-ray emission is possible through several processes. Often radiative wind 
shocks are considered, yet with no emission below 2 keV detected from WR 142, a line-driven 
flow instability model predicting soft X-rays from radiative wind shocks that is typically 
assumed for O stars (Gayley & Owocki 1995) cannot explain the hard emission. We discuss 
wind kinetic energy conversion, colliding wind shocks, magnetically confined wind shocks 
(MOWS), and coronal emission as possible thermal emission processes. 

5.4- 1- Kinetic Energy of the Wind 

Due to the extreme mass loss rate, the kinetic energy in the wind of WR 142 is enough 
to account for the X-ray luminosity. The kinetic energy of the wind supplies log L^ind ~ 
38.2 (ergs s~^) which is much greater than the X-ray luminosity observed of log Lx ~ 31 
(ergs s~^). The effects of clumping in the wind on the mass loss rate would not change this 
conclusion, reducing log L^ind by 0.6 dex (ergs s~^) at most. Thus even if the efficiency for 
converting the wind kinetic energy into X- rays were small, the wind still could supply the 
energy. 

5.4-2. Colliding Wind Shocks 

Because WR 142 has no known companion, a coUiding wind shock model does not 
obviously apply. However, such models could explain the high temperatures of the X-ray 
emission produced from the high wind speed of WR 142. Assuming v^ = 5500 km/s 

(Kingsburgh et al. 1995), the maximum X-ray temperature from a colliding wind shock (eq. 
1 of Luo et al. 1990) would be kT^ax ~ 1-95 keV ^ [v±/ (1000 km/s'^)]'^ > 59 keV for WR 
142. The symbol /j, is the mean atomic (ionic) weight per particle (ions -|- electrons). This 
approximation does not assume specific WO-type wind abundances (though // > 1 for WR 
stars), but shows that a colliding wind shock onto an unseen companion or its wind could 
explain the high temperatures indicated by APEC spectral fits. 

If we assume a BOV companion with a radius R* = 7.4 R0 (Allen 1976), we find a 
separation of ~ 1 AU is required to produce the observed Chandra X-ray luminosity through 
coUiding wind shocks (eq. 81 of Usov 1992). If we account for clumped winds, the separation 
scales as M^/^ and would reduce to 0.7 AU if the mass loss rate is reduced by a factor of 
2. At a distance of 1.23 kpc, a 1 AU separation between WR 142 and its hypothetical 
companion would produce an angular separation of much less than 1" (ft! 1 mas) and would 
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not be resolvable by Chandra. However, this separation is smaller than the escape radius 
for X-rays, which is 375 - 1500 R© (Section 5.1). A favorable orbital phase and geometry 
with the wind blown cavity around the companion along our line of sight would allow for 
the X-rays to escape (e.g. as occurs in 7^ Vel; Schild et al. 2004), especially with clumped 
winds. To consider shock emission from the fast winds impacting circumstellar material, 
more information on the CS material's density, geometry, and distance from WR 142 would 
be necessary. 



5.4-3. Magnetically Confined Wind Shocks 

For luminous stars, a surface magnetic field may be able to channel wind flow into shock 
collisions with sufficient velocity to produce hard X-ray emission. Generally, the magnetic 
field lines are thought to bend radiatively driven winds towards the magnetic equator, where 
the hemispheric winds collide. Such magnetically-confined wind theories have been used to 
explain X-ray emission for magnetic Ap-Bp stars (Babel & Montmerle 1997a) and young 
0-type stars (Babel & Montemerle 1997b, 9 Ori C in Gagne et al. 2005). The winds were 
not radiatively driven in the Ap-Bp model but the resulting shocks are unaffected. We can 
determine the minimum magnetic field strength necessary to confine the winds in comparison 
to the wind magnetic confinement parameter using equation 7 of ud-Doula & Owocki (2002), 
which is dependent upon the mass loss rate, Vqo, and radius of the star. We approximate 
the radius of WR 142 to be one solar radius (Abbott et al. 1986). We find that B ^/r^ 
20 kG at the stellar surface. Here, rj^, is the wind magnetic confinement parameter and is 
1 for marginal confinement, 10 for strong confinement that produces shocks strong enough 
to produce 1 keV emission (ud-Doula & Owocki 2002). Wind clumping could reduce the 
minimum magnetic field strength by a factor of ^/2 to 2. It remains to be shown that WR 
stars have such strong magnetic fields. 

Additionally, it can be shown that the maximum temperature at the shock front from 
MGWS could reach maximum temperatures of tens of keV (eq. 4 of Babel & Montmerle 
1997a) for high-velocity ionized metal-rich winds (// > 1.33) such as that of WR 142. The 

highest temperature is only possible at the shock front and will drop off to lower temperatures 
for most of the shocked plasma, but still easily could explain the observed hard emission from 
WR 142. It is important to note though that MGWS theories have not yet been extended to 
WR stars, where several problems such as larger mass loss rates have not yet been addressed. 
Additionally, a magnetic rotator model with modulated the winds was invoked to explain 
the periodic variations of the X-ray emission from 9^ Orionis G (Babel & Montmerle 1997b) 
but WR 142 here shows no such periodic X-ray emission or B-field detection so far. 
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5.4-4- Coronal Emission 

As suggested for some OB stars, X-rays could be formed in a thin corona at the base 
of the wind (e.g. Cassinelh & Olson 1979). However, our calculations suggest that the kT 
~ 2 keV X-rays emerge from very large radii, roughly r > 1000 R^., assuming a spherical 
homogenous wind. This is too far out to be attributed to a corona. Also, X-ray variability 
would be expected with coronal emission as magnetic (coronal) emission is usually associated 
with short-term X-ray variability (e.g. flares). No significant variability or flaring was 
detected in the Chandra observation. 



5.5. Nonthermal Emission Processes 

Electrons in the stellar wind accelerated to relativistic energies by the Fermi mech- 
anism and trapped in weak magnetic fields could scatter stellar UV photons to produce 
X-ray emission and 7-ray emission (Chen & White 1991a, b), as suggested for WR stars 
by Pollock (1987). The accelerated electrons undergo adiabatic cooling between shocks and 
thus can only lose heat through nonthermal processes, such as inverse Compton, nonther- 
mal bremsstrahlung, or synchrotron radio emission (Chen & White 1991a). Such processes 
may be at work in WR 142: X-ray spectra of WR 142 can be fitted with a power-law and 
Berkeley 87 lies close to an extended region of high-energy 7-ray emission in Cygnus (Abdo 
et al. 2007, 2009; Smith et al. 2005), though the primary source of the 7-ray emission has 
not yet been proven to be Berkeley 87 itself. 

For OB stars, the radiation will produce a power-law electron spectrum where the X- 
ray flux from inverse Compton emission should scale with energy as Fx oc E^^/^ for electron 
index n = 2 (Chen & White 1991a). Spectral fits in Table 3, when the given photon index 
Tpow = 1 is converted to an energy index, give Fx oc E°, though the photon index was poorly 
constrained. If the photon index is frozen to F = 1.5 so as to produce Chen & White's energy 
power-law index, the fit is reasonable but exhibits slightly more absorption and higher X-ray 
fiux and luminosity (Table 3 footnotes). 

The X-ray emission is predicted to be hard without large scale variations, and if present, 
variations only on timescales of hours to days (Chen & White 1991a). However, the Chen 
&; White (1991a) model is tuned to OB stars, such that inverse Compton is the assumed 
emission mechanism whereas bremsstrahlung may be just as important or even dominate 
inverse Compton for WR stars. Additionally, other factors, such as the coulomb effects in the 
dense WR wind, have not been fully evaluated (Chen & White 1991a). Still, the nonthermal 
emission process is plausible for some WR stars, especially where nonthermal radio emission 
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is detected (Pollock 1987). WR 142, is as of yet, undetected in the radio. Cappa et al. 
(2004) determined an upper limit of <0.90 mJy for the 3.6 cm flux density from WR 142. 
Despite a lack of detection of nonthermal radio emission at 3.6 cm, radio observations at 
longer wavelengths where nonthermal emission could lead to higher flux densities would be 
useful. The hard and essentially featureless X-ray spectrum of WR 142 could be nonthermal, 
but a higher signal-to-noise X-ray spectrum will be needed to distinguish between thermal 
and nonthermal emission. 



6. Comments on Wolf-Rayet Steir X-ray Emission 

The X-ray detection of WR 142 presents a major question: why don't we observe X-rays 
from the WC subtype when wc detect hard emission from a WO star? When plotting Lboi 
versus Lx, WR 142 lies in the region occupied by undetected WC stars (Figure 5). WOs are 
more evolved than WCs, and WO winds are even more metal-rich and should thus be more 
efficient absorbers of soft X-rays. 

However, WR 142 has a very high terminal wind speed and a very high effective temper- 
ature (Teff) compared to WC stars (Crowther 2007). A higher effective temperature would 
lead to a more fully ionized wind, which lowers the wind opacity {k) to X-rays (see eq. 6 of 
Ignace et al. 2000). The higher wind speed leads to a smaller radius of optical depth unity 
for X-rays of a given energy where R(t = 1, E) = M/(47rt;oo) x k (E) (eq. 9 of Ignace et 
al. 2000). In the case for WR 142, the terminal wind speed is greater and the wind opacity 
may be smaller than for WC stars, leading to a smaller R(r = 1, E). Thus, X-rays may be 
able to escape from smaller radii closer to the star, which could explain why WR 142 has 
been detected in the X-rays but why WC stars have (so far) not been detected. 

Another possible explanation is that WR 142 is not a single star (see Section 5.2). 
Some WC stars that are members of binary systems have been known to display hard X- 
ray emission (e.g. 7^ Vel, WC8+07.5) but the X-ray emission from WR 142 is much less 
luminous. WR+0 binaries display much stronger emission. Perhaps X-ray emission from a 
binary decreases as a WR star evolves from WN to WC to WO because of changes within the 
star or its wind. Maybe the luminosity difference lies in evolution of the unseen companion 
star. WR 142, if a binary, may not have an O star companion but some other less massive 
non- degenerate stellar companion. As discussed in Section 5.4.2, a close companion at the 
right separation could account for a lower X-ray luminosity using a colliding wind model. 

The hard X-ray emission and fast wind from WR 142 may indicate a colliding wind shock 
that could be explained by an as of yet undetected companion at close separation, such as a 
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BOV star companion at ~ 1 AU from WR 142 (Section 5.4.2). This close separation would 
not be resolved by Chandra. If a less massive unseen companion is present, then Chandra 
may be observing colliding wind shock emission at or near the surface of the companion 
(Usov 1992) and no intrinsic emission from the stars themselves. Even if the undetected 
companion star had no wind itself, the overwhelming wind from WR 142 would shock onto 
the companion's surface and produce the same effect. For example, a colliding wind shock 
onto a less massive companion could occur if the companion were an X-ray faint B- or A-star 
and any intrinsic X-ray emission from WR 142, if any exists, were completely absorbed, e.g. 
by its wind. If this is the case, then we do not detect X-ray emission from WR 142 itself 
and there is no contradiction with non-detections of single WC stars. 

Another thought to consider could involve the very high wind speeds of WR 142. What 
if the winds give rise to some exotic nonthermal emission processes, such as proposed by Chen 
& White (1991a), that do not get triggered in WR stars with lower terminal wind speeds or 
lower Teff? The detection of X-rays from the WO type WR 142 may even suggest that the 
absence of X-rays from WC type stars be considered tentative until sensitive observations of 
a larger sample of WC stars are obtained. 

7. Conclusions 

1. Chandra has detected hard, heavily- absorbed X-ray emission from the rare WO- type 
star WR 142. No soft emission below 2 keV was detected by Chandra. 

2. The observed X-ray flux is consistent with the flux from a previous XMM-Newton 
observation, within the uncertainties, and no significant X-ray variability was detected 
during the Chandra observation. 

3. Due to the faint emission from WR 142, lack of prominent emission lines, and low 
numbers of counts, statistics are unable to distinguish between thermal and nonthermal 
X-ray emission mechanisms when fitting the spectrum and both have been considered. 
If the emission is thermal, very high plasma temperatures are inferred. 

4. In addition to WR 142, Chandra detected seven luminous OB stars in Berkeley 87, 
while eight other B stars were undetected. The hard X-rays and excess absorption 
of WR 142 contrast with two X-ray bright OB stars in Berkeley 87, which display 
predominately soft X-ray emission and absorptions consistent with optical values. 

5. Though the X-ray emission mechanism in WR 142 is unclear, the hard X-ray spectrum 
observed by Chandra could be explained by a colliding wind shock onto an as yet 
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unseen companion, though the escape of the X-rays may require a geometry with the 
companion in front of WR 142. A coUiding wind shock interpretation would also resolve 
an apparent contradiction of non-detections in the X-rays of single WC stars while WR 
142 (a WO star) was detected. A key point not to be overlooked is that colliding winds 
can produce very hot plasma (as observed) even in the absence of magnetic fields. 

6. Alternatively, mechanisms that assume stellar magnetic fields such as MOWS or inverse 
Compton scattering (as formulated by Chen & White 1991a) could play a role in the 
X-ray emission of WR 142. However for MCWS, very strong surface fields of tens of 
kG would be required to effectively confine the powerful WR 142 wind (Sec. 5.4.3), 
also noted by Oskinova ct al. (2009) in their MCWS interpretation. There is so far 
no observational evidence of such strong B-ficlds in WR stars, nor do the existing X- 
ray data show any obvious signatures of impulsive X-ray flares that often accompany 
magnetic reconnection. If WR 142 has even a weak surface magnetic held, then X- 
ray production via inverse Compton scattering provides an attractive alternative to the 
extreme conditions required for magnetic wind conflnement. A key question is whether 
the X-ray spectrum of WR 142 is indeed a power-law, as expected for inverse Compton 
scattering. To answer this question, a higher signal-noise X-ray spectrum from a much 
deeper observation will be required. 
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Table 1. Stellar Properties of WR 142 



Parameter 


Value 


Notes^ 


Spectral Type 


W02 


(1)(2) 


Age (Myr) 


~3 - 5 


(3) 


d (pc) 


1230±40 


(4) 


V (mag.) 


13.4 


(2) 


Av (mag) 


6.1 


(5) 


J,H,Ks (mag) 


9.54,8.89,8.60 


(6) 


log M/Mq 


0.9±0.2 


(7) 


M (Mq yr-l) 


1.7 X 10-5 


(8) 


Voo (km s~^) 


5500±200 


(2) 


log L„ind (ergs s-l) 


38.2 


(9) 


log Lboi/L© 


5.35, 5.65 


(10) 



'^Refs. and Notes; (1) Crowther et al. 
(1998) (2) Kingsburgh et al. (1995) (3) Massey 
et al. (2001), Polcaro & Norci (1998), Turner ct al. 
(2010) (4) Turner et al. (2006) (5) van der Hucht 
2001, using Ay =0.9A„ (6) 2MASS Point Source 
Catalog (7) Smith et al. (1994) (8) Unclumped 
value from Barlow (1991), Clumping could reduce 
M by a factor of as 2 - 4 (Crowther 2007) (9) Pol- 
caro et al. (1991), also L^ind =(l/2)Mv^ (10) 
Determined by fitting photometric observations 
using a range of clumped mass loss rates, Oski- 
nova et al. 2009 
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Table 2. 


X-ray Properties 


of WR 142 and OB stars'^ 


in Berkeley 87 




Name 


R.A. 


Decl. 


Net Counts 




<B> 


Pconst 


Sp. Type 


Identification (offset) 




(,12000) 


(J2000) 


(cts) 


(koV) 


(koV) 


KS/GL 




(arcscc) 




HD 229059 


20 21 15.37 


-1-37 24 31.2 


341±19 


1.09,1.36,1.67 


1.54 


0.40/0.93 


Bl la*" 


HST 202115.36-1-372431.1 


(0.16)<^ 


Berkeley 87-4 


20 21 19.26 


-H37 23 24.6 


63±8 


0.96,1.19,1.67 


1.53 


0.10/0.49 


BO. 2 III 


HST 202119.25-H372324.7 


(0.16) 


Berkeley 87-13 


20 21 31.60 


-H37 20 40.6 


5±2 


...,1.25,... 


1.60 


0.66/... 


BO. 6 III 


HST 202131. 55-H372041. 4 


(1.00) 


Berkeley 87-24 


20 21 38.21 


-H37 25 17.4 


9±3 


1.69,1.99,3.49 


2.41 


0.49/... 


Bl lb 


HST 202138.33-H372616.3 


(1.80) 


BD-l-36 4032 


20 21 38.69 


-1-37 25 15.4 


254±16 


1.08,1.34,1.64 


1.43 


0.01/0.62 


08.5 III'' 


HST 202138.68-1-372616.2 


(0.23) 


Berkeley 87-26 


20 21 39.72 


-1-37 25 05.6 


11±3 


0.92,1.09,1.78 


2.00 


0.70/... 


BO. 5 lab 


HST 202139.75-1-372505.9 


(0.47) 


WR 142 


20 21 44.35 


-1-37 22 30.7 


46±7 


2.99,4.28,5.32 


4.31 


0.77/0.40 


W02'' 


HST 202144.36-1-372230.5 


(0.23) 


Berkeley 87-32 


20 21 47.35 


-H37 26 32.2 


11±3 


0.95,1.09,1.57 


1.33 


0.91/... 


BO. 5 III 


HST 202147.35-H372632.0 


(0.18) 



Note. — X-ray data using events in the 0.3 - 8 keV range within a 3<7 wavdetect ellipse, which for WR 142 is 11.25 pix^. Tabulated quantities are: 
J2000.0 X-ray position (R.A., decl.); net counts and net counts error from wavdetect (accumulated in a 70148 s exposure, rounded to the nearest integer, 
background subtracted and PSF-corrected) ; 25%, 50%, and 75% photon quartile energies (E25, E50, and E75), mean photon energy (<E>); probability of 
constant count-rate determined by the Kolmogorov-Smirnov (KS) test and the Gregory-Loredo (GL) algorithm, where Pconst ^ 0.05 would indicate likely 
variability; spectral type as indicated by SIMBAD or otherwise noted; and HST GSC candidate counterpart identification within a 2" search radius. The 
offset (in parentheses) is given in arc seconds between the X-ray and counterpart position. Ellipses indicate insufficient counts for reliable measurement. 

^Faint B stars not detected: Berkeley 87 nos. 7 (B), 9 (B0.5 V), 15 (V439 Cyg, BOe/B[e] - references in Sec. 4.3), 16 (B2 V), 18 (Bl V), 31 (Bl V), 34 
(B), and 38 (B2 111). Spectral types are from the SIMBAD database. 

'^Spectral type from Massey et al. (2001). 

"^Two close HST sources have nearly the same sexagesimal coordinates and are indistinguishable by Chandra. 
"^Spectral type references in Table 1. 
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Table 3. Chandra Spectral Fits for WR 142 



Parameter 



Model 


IT APEC 


POWER 


POWER 4 


Nh (1022 cm-2) 


5.66 


4.61 


4.55 


kT (keV) 


{20.}- 






norm (lO"") 


6.36 






r b 

J- pow 




1.0 


{1.0} 


normpow (10^*') 




2.31 


2.32 


Elinc (kcV) 






{6.67} 


0-line (keV) 






{0.05} 


normiine (10"^ ) 






8.68 


Abundances'^ 


WO 






C-statistic'' 


215.90 


215.08 


214.73 


Fx (lO-l-* ergs cm-2 s"!) 


1.65 (3.66) 


1.74 (2.93) 


1.77 (2.96) 


Fx,iine (10-" ergs cm'^ s"!) 






8.66 (9.28) 


log Lx (ergs s-^) 


30.82 


30.72 


30. 73 



Note. — Based on fits of ACIS spectra using XSPEC vl2.5.0. Tlic data were un- 
binned and not background subtracted (background is negligible). The tabulated pa- 
rameters are absorption column density (Nh), plasma temperature (kT), XSPEC nor- 
malization (norm), photon power-law index (Fpow), power-law model normalization 
(normpow), Gaussian line centroid energy (Eiino), line width {un^^ = FWHM/2.35), 
line normalization (normiino). Quantities enclosed in curly braces were held fixed dur- 
ing fitting. WO abundances arc from Tabic 1 of van dcr Hucfit ct al. (1986), except 
the hydrogen and nitrogen abundances arc arbitrarily set to the small non-zero value 
1 X 10"® for compatibility with XSPEC. Solar values from Anders & Grevesse (1989) 
are used for any elements in the XSPEC abundance table that are not listed in van 
der Hucht et al. (1986). X-ray flux (Fx) is the observed (absorbed) value followed 
in parentheses by the unabsorbed value in the 0.3 - 8.0 keV range. The continuum- 
subtracted Gaussian line flux (Fx.unc) is measured in the 6.5 - 6.9 keV range. X-ray 
luminosity (Lx) is the unabsorbed value in the 0.3 - 8.0 keV range. A distance of 
1230 pc is assumed (Turner et al. 2006). 

'^If kT is allowed to vary, it runs away to the maximum temperature of 64. keV 
allowed by the APEC model. If frozen at a higher temperature of kT = 50 keV, the 
fit is similar and converges to Nh = 5.38 X 10^^ cra~^, C-statistic = 215.56, Fx = 
1.67 (3.43) X 10-" ergs cm-^ s"!, and log Lx = 30.79 ergs s"!. 

''A photon index Fpow equates to Fx oc E-^^p""-^). That is, values of Fpow > 1 
correspond to flux decreasing with increasing energy E. 

•^If the photon index is frozen to Fpow = 1-5, the fit gives Nh = 5.91 X 10^^ cm^^, 
C-statistic = 215.73, Fx = 1.64 (4.03) x IQ-^* ergs cm-^ s'^, and log Lx = 30.86 
ergs s"^. 

'^Varying the abundances did little to improve fits. 

''Regrouped spectra to a minimum of 6 counts per bin enabled the use of statis- 
tics. For IT APEC models, this yielded a lower 90% confidence bound on kT>3.6 
keV. 90% bounds on Nh were in the range [2.8 - 9.6] X 10^^ cm-^ for models frozen 
at kT = 20 keV. We were not able to find stable confidence bounds on Fpow for the 
power-law models. 
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Table 4. Chandra Spectral Fits for X-ray Bright OB Stars 



Parameter 


BD+36 4032 


HD 229059 


Berkeley 87-4 


Model 


IT APEC 


IT APEC 




IT APEC 


Abundances 


solar 


solar 




solar 


Nh (1022 cm-2) 


1.13 [0.99 - 1.45] 


1.16 [1.01 - 


1.24] 


0.82 [0.59 - 1.16] 


kT (keV) 


0.56 [0.35 - 0.69] 


0.64 [0.47 - 


0.81] 


0.54 [0.27 - 0.74] 


norm (10"*) 


1.16 [0.69 - 4.31] 


1.47 [0.76 - 


1.75] 


0.20 [... - 1.51] 


x2/dof 


23.7/18 


30.7/28 




2.8/4 




1.32 


1.10 




0.70 


t'x (lO^^** ergs cm~^ ^~^) 


0.20 (2.91) 


0.32 (3.93) 




0.05 (0.53) 


log Lx (ergs s~^) 


31.72 


31.85 




30.98 



Note. — Based on fits of ACIS spectra binned to a minimum of 10 counts per bin using 
XSPEC vl2.5.0. Square brackets enclose 90% confidence intervals. An ellipsis means the 
algorithm used to compute 90% confidence limits did not converge. The X-ray flux is the 
absorbed value, followed in parentheses by the unabsorbed value. A distance of 1.23 kpc 
is assumed. Solar abundances are from Anders & Grevesse (1989). 
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WR 142 




□ - V439 Cyg 



HD 229059 



Berkeley 87-4 - 



22:30.0 -WR142 

45.0 - . 40.0 



Fig. 1. — Left: An ACIS-I image (0.3 - 8 keV, log intensity scale) centered near WR 142. 
Indicators give the WR 142 position as follows: circle marks HST GSC 2.3.2 (optical), plus 
marks USNO-Bl.O (optical), and the cross marks both the Chandra and 2MASS (near- 
infrared) positions. WR 142 was only 28" from the aimpoint. The 3a ellipse determined 
by wavdetect is shown and encloses 46 ± 7 net counts. Right: An AClS-1 wide field image 
(0.3 - 8 keV, log intensity scale) showing the X-ray bright OB stars relative to WR 142, 
where circles mark the Chandra positions and the square marks a 2MASS position for the 
non-detection of V439 Cyg. All coordinates are J2000.0 and pixels are 0.492". 
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Fig. 2. — A histogram showing the observed X-ray event energies for WR 142 in the energy 
range from 0.3 - 12 keV {wavdetect Sa eUipse extraction region). The bin size is 1 keV. 



o 
o 



- BD+36 4032 (08.5III) 




0.5 1 2 5 

Energy (keV) 

o c — ' ' ' ' ' ' — I : ' ' ' 

o : . 



- HD 229059 (B1 la) 



1 

o - 



WR 142 (W02) 




o r 



0.5 1 2 5 

Energy (keV) 

Fig. 3.— Chandra ACIS-I spectra of WR 142 (46 counts) and OB stars, top: the 08.5 III 
star BD+36 4032 (249 counts), bottom: the Bl la star HD 229059 (341 counts). Spectra 
were rebinned to a minimum of 5 counts per bin. The straight hues between data points are a 
histogram and not a fitted model. The spectra show a clear contrast. The WR 142 spectrum 
is dominated by hard emission above 2 keV while soft emission below 2 keV dominates the 
OB stars. 
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Fig. 4. — A Chandra lightcurve for WR 142 in the 0.3 - 8 keV energy range {la error bars). 
The extraction region was a 3a eUipse (Figure 1) given by wavdetect output. 
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Fig. 5. — Adapted from Skinner et al. (2010b), with new points for WR 142 and WR 16. 
Note that WR 142 lies in the region of undetected WC stars. The dashed line marked 10~^ 
approximately follows the Lx versus Lboi relation for O- and early B-stars from ROSAT 
(Berghofer et al. 1997), though a large scatter exists. The Lx value for WR 16 is based on 
a recent XMM-Newton observation (Observation Id 0602020301, obtained on 2009 Dec. 28). 
EPIC MOS images show a faint excess above background at the WR 16 position (13 net 
M0S2 counts in the 0.3 - 8 keV range, r = 15" extraction circle = 68% encircled energy, 33 
ksec of usable M0S2 exposure, 0.394 counts ksec^^). Using PIMMS, the M0S2 count rate 
equates to an unabsorbed X-ray luminosity log Lx(0.3 - 8 keV) = 30.82 ergs s~^, assuming 
d = 2.73 kpc and Ay = 1.85 (van der Hucht 2001), and an intrinsic 2T Raymond-Smith 
thermal plasma source spectrum with kTi = 0.6 keV and kT2 = 3.5 keV. Because the X-ray 
detection of WR 16 is of low significance, its Lx value is plotted as an upper limit. Lboi for 
WR 16 is from Hamann et al. (2006). 



